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FOREWORD 
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14670246,  "Automated  Design  of  Advanced  Aerospace  Structures".  The 
work  was  carried  out  in  the  Design  and  Analysis  Methods  Group  of  the 
Analysis  and  Optimization  Branch  (FBR),  Structural  Mechanics  Division, 
Air  Force  Flight  Dynamics  Laboratory  (AFFDL),  Wright-Patterson  AFB, 
Ohio. 


The  time  period  of  the  effort  was  June  1975  - June  1976. 


The  manuscript  was  released  in  August  1976. 


AFPDL-TR- 76- 1 1 7 


TABLE  OF  CONTENTS 


SECTION 

PAGE 

I 

INTRODUCTION 

1 

II 

OPTIMALITY  CRITERIA  AND  RECURRENCE  RE! AT  ION 

3 

III 

DESIGN  PROCEDURE 

8 

IV 

ILLUSTRATIVE  EXAMPLE 

11 

1.  Stress  Constraint  Problem  (Conventional  Design) 

12 

2.  Twist  Constraint  Problems 

12 

V 

CONCLUSION 

14 

REFERENCES 

39 

v 


PRECEDING  PAGE  BLANK- NOT  FILMED 


r 


AFFDL-TR-76-117 


ILLUSTRATIONS 


FIGURE 

PAGE 

1 

Planform  of  Wing 

15 

2 

Definition  of  Twist 

15 

3 

Type  of  Elements 

16 

4 

Model  of  Wing 

17 

5 

Optimum  Design  of  Wing  (Stress  Constraint) 

18 

6 

Iteration  for  Twist  Constraint  (Wash  out  -7°) 

19 

7 

Optimum  Design  of  Wing  (Twist  Constraint,  Wash  out  -7") 

20 

8 

Iteration  for  Twist  Constraint  (Wash  in  2°) 

21 

9 

Optimum  Design  of  Wing  (Twist  constraint,  Wash  in  2°) 

22 

vi 


W 


AFFDL-TR- 76- 1 1 7 


TABLES 


TABLE 

PAGE 

1 

Element  Connections 

23 

2 

Coordinates  of  Nodes  and  Applied  Loads 

24 

3 

Properties  and  Allowable  Stresses 

26 

4 

Iterative  History  for  the  Stress  Constraint  Problem 

27 

5 

Element  and  Layer  Thicknesses  for  the  Minimum 
Weight  Stress  Constraint  Design 

28 

6 

Shear  Panel  and  Bar  Dimensions  for  the  Minimum 
Weight  Stress  Constraint  Design 

30 

7 

Deflections  at  Tip  for  the  Minimum  Weight  Stress 
Constraint  Design 

30 

8 

Iterative  History  for  the  Twist  Constraint 
Problem  (Wash  out  -7°) 

31 

9 

Deflections  at  the  Tip  for  the  Minimum  Weight 
Twist  Constraint  Design  (Wash  out  -7°) 

31 

10 

Element  and  Layer  Thicknesses  for  the  Minimum 
Weight  Twist  Constraint  Design  (Wash  out  -7°) 

32 

11 

Shear  Panel  and  Bar  Dimensions  for  the  Minimum 
Weight  Twist  Constraint  Design  (Wash  out  -7°) 

34 

12 

Iterative  History  for  the  Twist  Constraint 
Problem  (Wash  in  2°) 

35 

13 

Deflections  at  Tip  for  the  Minimum  Weight  Twist 
Constraint  Design  (Wash  in  2°) 

35 

14 

Element  and  Layer  Thicknesses  for  the  Minimum 
Weight  Twist  Constraint  Design  (Wash  in  2°) 

36 

15 

Shear  Panel  and  Bar  Dimensions  for  the  Minimum 
Weight  Twist  Constraint  Design  (Wash  in  2°) 

38 

L 


AFFDL-TR-76-117 


SECTION  I 
INTRODUCTION 

During  the  last  few  years  the  design  concept  of  aeroelastic  tailoring 
aircraft  wing  structures  has  gained  considerable  interest.  The  basic 
idea  of  this  concept  is  to  control  the  anisotropic  properties  of  the  wing 
structure  so  that  the  wing  will  have  deformation  characteristics  under 
load  such  as  twist  and  camber  which  provide  beneficial  aerodynamic  per- 
formance while  also  satisfying  strength  and  flutter  requirements.  The 
elastic  properties  of  the  composite  laminate  can  be  defined  by  a suitable 
selection  of  the  number  of  layers,  the  fiber  direction  in  each  layer  and 
the  number  of  plies  in  the  layer  having  the  same  fiber  orientations.  In 
this  report,  it  is  assumed  that  the  twist  of  the  lifting  surface  is  of 
primary  importance  and  the  structure  is  to  be  designed  for  a specific 
value  of  the  twist  between  two  points  on  a chord  and  satisfy  the 
strength  requirements. 

Methods  developed  on  the  basis  of  optimality  criteria  are  found  to 
be  efficient  and  have  received  wide  attention  in  recent  years.  An  over- 
all review  of  these  methods  is  given  in  References  1 and  2.  The  objective 
in  optimization  of  layered  composites  treated  in  this  report  is  to  find 
the  thickness  of  each  layer  in  the  composite  elements  such  that  the 
designed  structure  satisfies  the  strength  criteria  in  all  the  elements 
and  achieves  the  specified  displacement  pattern  as  well.  The  problem 
of  designing  a wing  type  structure  with  specified  twist  and  displacements 
is  discussed  in  Reference  3.  In  this  reference,  the  twist  is  specified 
in  terms  of  two  unequal  displacement  constraints  at  two  points  on  a i herd 
in  the  structure  and  the  problem  is  treated  as  a multiple  displacement 
constraint  problem.  In  the  method  of  Reference  3,  the  twist  constraint 
could  not  be  separated  from  the  displacement  constraints.  In  the 
approach  discussed  in  this  report,  however,  the  restriction  of  specifying 
the  displacements  is  removed  and  one  can  specify  tin’  required  rotation 
of  a line  between  any  two  points  on  a chord. 

The  planform  of  a wing  structure  is  shown  in  I iguro  1 . A and  It  are 
points  on  the  trailing  and  leading  edge  respectively,  at  the  tip  of  the 
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wing.  In  Figure  2,  AB  is  the  undeflected  position  of  the  tip  of  the 
wing.  In  Figure  2(a),  is  the  deflected  position  of  the  tip  for 

the  conventional  design.  The  conventional  design  is  defined  as  a design 
that  satisfies  the  stress  constraints  only  and  has  no  specific  deflection 
or  twist  requirement.  In  the  conventional  design  the  apparent  center  of 
rotation  is  in  front  of  the  leading  edge  and  the  tip  rotates  through 
G1  ('wash  out'  condition).  Two  twist  requirements  may  be  defined  either 
(1)  to  increase  the  conventional  design's  wash  out  twist  as  shown  in 
Figure  2(b)  so  that  ©2  > ©i  or  (2)  to  reverse  the  direction  of  the 
conventional  design's  twist  as  shown  in  Figure  2(c).  The  second  require- 
ment corresponds  to  the  'wash  -in'  condition.  In  the  first  case,  the 
apparent  center  of  rotation  moves  towards  the  wing  point  B,  and  in  the 
second  case  the  center  of  rotation  moves  in  back  of  the  trailing  edge. 

In  subsequent  sections  for  the  convenience  of  discussion,  the  'wash  out' 
and  'wash  in'  conditions  will  be  referred  to  as  negative  twist  and 
positive  twist  respectively. 


In  this  report,  the  application  of  the  optimality  criteria  method 
is  illustrated  for  the  above  twist  conditions  by  designing  a wing  structure. 
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SECTION  II 

OPT I MI L I T Y CRITERIA  AND  RECURRENCE  RELATION 


The  structure  is  discretized  into  m finite  elements  where  each 
element  contains  n subelements  corresponding  to  the  number  of  layers  and 
each  layer  consists  of  a number  of  plies  with  the  same  fiber  orientation. 
The  relation  between  the  applied  load  vector  { R}  and  the  displacement 
vector  { r}  can  be  written  as 

f K ] { r } = {R}  (1) 


where  [K]  is  the  stiffness  matrix  of  the  total  structure.  The  nodal 
displacements  f r}^  of  the  ith  element  and  the  nodal  forces  { f}.  . of  the 
jth  layer  of  the  ith  element  are  related  by 


Ik].,  {r}  = {f}.. 

1 1J 


(2) 


where  [k]^j  is  the  stiffness  matrix  of  the  jth  layer  of  the  ith  element. 
The  ith  element  strains  ( el.  are  related  to  the  nodal  displacements  by 


{c)i  = [c]  {r}. 


(3) 


where  [c]  represents  the  matrix  of  strains  in  the  ith  element  for  a 
unit  displacement  vector.  The  stress  strain  relation  for  the  jth  layer 
of  the  ith  element  is 


[o],.  = [ D ] . .{el. 

1J  1 


(A) 


where  [D]  is  the  matrix  of  elastic  constants  and  [o]^  are  the  stresses 
in  the  jth  layer  of  the  ith  element. 


The  total  weight  of  the  structure  is  given  by 

m n 

W = 2 fi 


1=1 


j = l 


1 . . A,  . 
iJ  ij 


(5) 


where  . is  the  density,  A^j  is  the  thickness  of  the  jth  layer  of  the 
ith  element,  and  is  the  area  of  the  ith  element. 
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Let  the  prescribed  twist,  between  two  points  in  a structure  be 
denoted  by  c.  Then  the  constraint  equation  can  be  written  as 

(ri*  (Si  - c = 0 


(6) 


where  r^  represents  the  transpose  of  the  displacement  vector  and  rS 
represents  the  virtual  load  vector.  The  latter  vector  consists  of  forces 
corresponding  to  a unit  couple  applied  between  the  two  points  where  the 
twist  is  specified.  Using  Equations  5 and  6,  the  Lagrangian  F can  be 
written  as  _ m n 

p = j Aij  + *:s}  - Cl  = 0 (7) 

r-J  i = l -> 


where  ■ is  the  Lagrange  multiplier.  Differentiating  Equation  7 with 


respect  to 

the 

design 

var 

iabl es 

Aij 

gives  the  necessary  conditions 

for 

an  optimum. 

p 

. . J.  . + 
ij  i 

A 

3 j 

3A.y  : 

r L 

{Si 

= 

0 

(8) 

. 

or 

qj 

- 

1 

1 

p , 

. . p . + 
ij  j. 

\ 

■ a'  . 

{ r : c 

f K 1 (s 

= 0 

(9) 

\ 

L tj 

- 

1 

j 

- - t . 2 , - - - , m 
j = 1,2,  - - n 

where  s is  the  displacement  vector  correspond! ng  to  the  virtual  load 
vector  Si. 


Differentiating  Equation  1 with  respect  to  the  design  variable  A.,  gives 

(10) 


>k  , rt  , a 

La  ,r  + lK')  , • fr}  = o 


1 I 


l 1 


Multiplying  Fguation  10  by  [K]  and  rearranging  gives 


r I ~ - 1 r • 


i l 


>K 

) A . . 
i.)  _ 


f K 1 


(11) 


Substituting  Equation  11  in  Equation  9 gives 


::  - - ' ^ ' J-  ....  0 

i I 


I 1 1 


(12) 


')  = 1.2, n 
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For  the  case  of  a constan 
to  A-  ■ , thus 

' J 


t strain  membrane  element, 


3K 

3A.  . 
ij 


1_ 

A.  . 
ij 


[k]jj  is  proportional 
(13) 


Hence,  Equation  12  can  be  written  as 


{ r 1 . [kl..  {s}. 

h - - 0 

1J 


(14) 


where f s}..  are  the  nodal  displacements  of  the  ith  element  due  to  the 
virtual  load  vector  {S}.  Hence,  the  optimality  condition  for  the  twist 
constraint  problem  is 


1 


{r fk]..  (si, 
> i. J-J 1 


(15) 


or 

e . . 

i = > 

0 . . 
■ IJ 


(16) 


1 * 1 , 2 , - - -■  m 
j * 1,2, n 


where 


i ) 


■ ; •• ',j 

A . . 1. 
i I i 


(I  /) 


e.  . is  the  virtual  strain  energy  density  of  the  jth  layer  of  the  ith 
element. 


The  * im.tl  i ty  condition  will  now  he  used  to  derive  a recurrence 
relation. 


The  design  variable  A. ^ can  be  written  as 

A . “ A» 

i i i j 

where  is  a constant,  scalinq  parameter  and  is  tiie  norma  I iced 
relative  design  vector.  In  the  case  of  an  isotrcgn  elo  . m , 1 and 
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represents  either  the  thickness  of  a plate  element  or  the  cross- 
sectional  area  of  a bar  element. 

Introducing  Equation  18  in  Equation  1 yields 

[K'l  { r ' } = { R } (19) 

where  [K1]  is  the  stiffness  matrix  for  the  total  structure  obtained  by 

using  the  normalized  relatives  sizes  a.,  of  the  elements  and 

i J 

{r ' } = A {r}  (20) 


At  the  element  level,  the  relations  between  the  actual  quantities 

and  the  relative  quantities  are 

{r ' } . = A lr 1 . 
i l 


{s’).  = A (s). 
l i 


(21) 


where  the  prime  denotes  the  relative  quantities.  Similarly,  the  relation 
between  the  relative  stresses  and  the  actual  stresses  in  an  element  can 
be  written  as 


{a'}..  = A (a) 

ij 


(22) 


where  the  relative  stresses  { o' 1.^  are  obtained  by  using  the  relative 
displacements  { r ' } ^ . 

From  Equation  22  it  is  seen  that  the  actual  stress  in  the  element  can 
be  modified  to  satisfy  the  selected  strength  criteria  by  selecting  a 
suitable  value  of  the  scaling  parameter  A. 


Introducing  Equation  21  in  Equation  16,  the  optimality  criteria 
can  be  written  as 


1 


itjj. 

A2  p 


1=1, m 

.1  = 1, n 


(23) 
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where  e'.,  is  obtained  from  Equation  17  by  replacing  all  the  actual 
quantities  by  their  relative  values.  Multiplying  Equation  23  by  A af  . 
and  taking  the  square  root  of  the  resulting  equation  yields 

vU2 


(a^A)  = B 


where  B is  a constant. 


(24) 


Equation ^24  can  be  written  in  an  iterative  form  as 


(a . . A)  = B (a . . ) 
lj  v+1  ij'vi 


Vl/2 


'ij 


(25) 


where  v+1  and  v refer  to  the  iteration  number.  Since,  the  design 
variables  are  expressed  as  normalized  relative  variables,  it  is  not 
necessary  to  evaluate  the  constant  B in  Equation  25  explicitly. 
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SECTION  III 
DESIGN  PROCEDURE 


A computer  program  based  on  the  finite  element  method  to  design  a 
minimum  weight  structure  with  a specified  twist  constraint  was  written 
incorporating  Equation  25.  Four  types  of  elements  are  included  in  the 
program:  (1)  constant  strain  triangles,  (2)  quadrilaterals  constructed 

from  four  constant  strain  triangles,  (3)  shear  panels,  and  (4)  bars. 

The  elements  are  defined  by  the  node  numbers  as  shown  in  Figure  3.  The 
triangles  and  quadrilaterals  may  be  layered  composite  elements  with 
specified  fiber  directions.  The  choice  of  the  number  of  layers  and  the 
fiber  orientation  is  arbitrary.  The  designer  can  select  any  appropriate 
number  of  layers  and  any  fiber  orientation  in  each  layer. 

The  basic  steps  involved  in  the  iterative  procedure  are:  (1)  analyze 

the  structure  using  the  relative  design  vector  and  find  the  nodal  displace- 
ments and  the  element  stresses,  (2)  find  the  scaling  parameters  which 
satisfy  the  strength  criteria  for  all  elements  and  the  twist  constraint, 

(3)  modify  the  relative  design  vector  by  using  Equation  25.  These 
basic  steps  are  repeated  for  a specified  number  of  iterations.  The  weight 
of  the  structure,  which  is  a good  indication  of  the  convergence  of  the 
iterative  procedure,  is  evaluated  during  each  iteration. 

The  two  scaling  parameters  are  As  for  satisfying  the  strenqth 
criteria  and  for  satisfying  the  twist  constraint.  The  design  with 
Aij  = As<i.j j will  satisfy  the  strenqth  criteria  for  all  the  elements. 

Since  the  scaling  parameter  is  the  same  for  all  the  elements,  for  some 
elements  the  strength  criteria  will  be  oversati sf i ed . The  design  with 
Aij  = A(jaij  satisfy  the  twist  constraint.  If  A^  >_  A$  then  the  design 
will  satisfy  the  twist  constraint  and  the  strenqth  criteria.  However,  if 
Ad  < A , then  the  design  with  the  scaling  parameter  A^  will  violate  the 
strength  criteria  in  certain  elements.  Since  the  objective  here  is  to 
satisfy  the  twist  constraint  and  the  strenqth  criteria,  the  acceptable 
designs  are  those  where  A^  > A . 


8 


AFFDL-TR-76-1 1 7 


A design  can  be  scaled  to  achieve  the  specified  twist  only  if  the 
structure  is  twisting  in  the  same  direction  as  the  direction  of  the 
specified  twist.  For  example,  if  the  unsealed  design  is  twisting  in  a 
negative  direction,  as  in  Figure  2(b),  then  the  design  cannot  be  scaled 
to  achieve  a positive  twist.  Similarly,  the  design  cannot  be  scaled  to 
achieve  a zero  twist,  (ie.  the  two  points  A and  B in  Figure  2 moving 
equal  distances  after  the  structure  deforms).  However,  a structure  can 
be  designed  for  zero  twist  by  interpolating  between  the  two  designs  which 
twist  in  the  positive  and  negative  directions  during  consecutive  iterations. 

The  twist  is  defined  as  the  angle  between  any  two  lines  each  con- 
necting the  same  two  points  in  a structure  before  and  after  deformation. 
Since  the  line  connecting  the  two  points  in  a structure  moves  in  a three 
dimensional  space,  the  correct  scaling  parameter  can  be  found  only  by 
numerical  iterations.  Direct  scaling,  as  defined  by  Equation  21  for 
displacements,  can  be  used  only  if  the  line  moves  parallel  to  one  of  the 
coordinate  planes  during  deformation. 

In  scaling  the  designs,  one  of  two  procedures  is  followed  depending 
on  the  direction  of  the  specified  twist  constraint.  If  the  wing  is  to 
be  designed  for  positive  twist,  i.e.,  wash  in  condition,  A^  is  first 
calculated  for  the  relative  design  vector  obtained  by  using  Equation  25. 
Then,  based  on  A^,  the  relative  sizes  of  the  passive  elements  are  changed 
so  that  the  stresses  in  these  elements  satisfy  the  strength  criteria 
identically.  The  passive  elements  include  those  with  negative  virtual 
strain  energy  in  the  last  iteration  and  also  those  which  violate  the 
strength  criteria  with  A^  as  the  scaling  parameter.  The  elements  with 
negative  virtual  strain  energy  cannot  be  modified  by  using  Equation  25. 

They  can  be  resized  by  selecting  their  relative  thickness  such  that  the 
stresses  in  these  elements  satisfy  the  strength  criteria  identically. 

After  this  modification,  the  structure  is  reanalyzed,  a new  A^  is  calcu- 
lated and  the  passive  members  are  remodified.  This  process  of  correcting 
Ad  and  modifying  the  passive  members  is  continued  for  a specified  number 
of  cycles  or  until  the  difference  between  the  previously  calculated  A^ 
and  the  new  A^  is  smaller  than  a prescribed  value.  In  the  present  compu- 
ter program  a maximum  of  three  cycles  is  allowed  for  the  adjustment  of 
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the  sizes  of  the  passive  elements.  After  three  cycles,  the  sizes  of  all 
elements  with  positive  virtual  strain  energy  are  changed  by  using 
Equation  25  and  the  sizes  of  all  elements  with  negative  virtual  strain 
energy  are  modified  based  on  their  stresses  from  the  last  iteration. 

In  the  second  procedure  where  the  wing  is  to  be  designed  for  negative 
twist,  i.e.  wash  out  condition,  A$  is  first  calculated.  In  determining 
this  value,  those  elements  with  negative  virtual  strain  energy,  the 
shear  panels  and  the  bars  are  not  taken  into  consideration.  Then  based 
on  As,  the  elements  with  negative  virtual  strain  energy  and  the  elements 
which  violate  the  strength  criteria  are  modified  so  that  they  satisfy 
the  strength  criteria  identically.  The  process  of  modifying  the  passive 
members  is  done  three  times.  The  scaling  parameter  is  then  calculated 
to  satisfy  the  specific  twist  constraint.  It  has  been  found  that  the 
number  of  cycles  required  to  modify  the  passive  members  decreases  as  the 
design  approaches  the  minimum  weight  design. 


in 
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SECTION  IV 

ILLUSTRATIVE  EXAMPLE 

In  order  to  illustrate  an  application  of  the  proposed  method,  a 
wing  structure  is  designed  for  (1)  stress  constraint,  (2)  negative  twist 
constraint  (wash  out  condition)  and  (3)  positive  twist  constraint  (wash 
in  condition).  The  structure  is  designed  with  stress  constraint  only  in 
order  to  establish  the  magnitude  and  direction  of  the  twist  of  the  tip 
of  the  wing  for  the  conventional  design. 

The  finite  element  model  of  the  wing  structure  is  shown  in  Figure  4. 
The  nodes  of  the  bottom  skin  have  odd  numbers.  The  structural  model 
has  88  nodes  and  158  members.  The  top  and  bottom  skins  are  idealized  by 
membrane  quadrilaterals  and  triangles.  The  spars  and  ribs  are  idealized 
by  shear  panels  and  posts.  The  element  numbers  and  the  node  numbers 
defining  the  elements  are  given  in  Table  1.  The  convention  used  to 
define  the  elements  is  given  in  Figure  3.  The  structure  is  subjected 
to  a single  loading  condition  as  given  in  Table  2.  The  coordinates  of 
the  node  numbers  are  also  given  in  this  table. 

The  skin  elements  consist  of  four  layers  with  fibers  in  the  0°, 

90°,  +45°  and  -45°  directions.  The  number  of  +45"  and  -45"  plies  may  not 
be  equal.  The  0°  fibers  are  parallel  to  the  direction  of  the  middle 
spar  as  defined  by  connecting  nodes  4 and  84  in  the  top  skin  or  3 and  83 
in  the  bottom  skin.  The  elastic  constants  and  the  allowable  strengths  for 
graphite  epoxy  and  aluminum  are  given  in  Table  3.  The  shear  panels  and 
posts  are  aluminum.  The  maximum  stress  criteria  are  used  to  size  the 
membrane  elements.  However,  only  the  stresses  in  the  fiber  direction 
and  the  shear  stresses  in  each  layer  are  checked  against  the  maximum 
allowables.  For  the  starting  design  the  relative  sizes  of  the  bar 
elements  are  1.0  and  the  thicknesses  of  the  plate  element  are  0.1.  The 
membrane  element  for  the  starting  design  are  assumed  to  have  an  equal 
percentage  of  plies  in  the  four  fiber  directions. 
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1.  STRESS  CONSTRAINT  PROBLEM  (CONVENTIONAL  DESIGN) 

The  structure  is  designed  by  using  the  recurrence  relation  given  in 
Reference  3.  This  recurrence  relation  is  similar  to  Equation  25,  except 
that  the  relative  virtual  strain  energy  term  e'.j  is  replaced  by  the 
relative  strain  energy  due  to  the  applied  set  of  loads. 

The  iteration  history  for  the  stress  constraint  problem  is  given  in 
Table  4.  The  initial  weight  of  the  structure  with  an  equal  percentage 
of  plies  in  the  four  fiber  directions  in  the  membrane  elements  is  312.84  lbs. 
The  least  weight  design  of  45.45  lbs  is  obtained  in  the  14th  iteration. 

The  sizes  of  the  elements  corresponding  to  the  minimum  weight  design 
are  given  in  Tables  5 and  6.  In  these  tables  the  quantities  represent 
the  thicknesses  for  the  plate  elements  and  the  cross-sectional  areas  for 
the  bar  elements.  Figure  5 shows  the  distribution  of  the  number  of 
plies  in  the  four  fiber  directions.  The  number  of  plies  is  obtained 
by  dividing  the  thickness  of  the  layer  by  the  thickness  of  a single 
ply  and  rounding  it  to  the  next  higher  number.  The  deflections  for  the 
final  stress  constraint  design  at  nodes  7 and  9 at  the  tip  of  the  wing 
are  given  in  Table  7.  The  tip  twists  through  an  angle  -5.458°.  This 
is  the  angle  between  the  lines  joining  nodes  9 and  7 before  and  after 
loading. 

2.  TWIST  CONSTRAINT  PROBLEMS 

The  wing  structure  is  designed  for  two  twist  constraint  conditions. 

The  purpose  of  the  two  examples  is  to  show  the  application  of  the  twist 
constraint  algorithm  (Equation  25)  in  designing  a structure  with  positive 
or  negative  twist.  The  two  angles  selected  for  the  twist  constraint 
problems  are  -7°  (wash  out)  and  +2°  (wash  in).  For  the  twist  constraint 
problem,  recall  that  there  are  two  scaling  parameters  A and  A^.  The 
designs  corresponding  to  A?  and  will  be  referred  to  as  'stress  design' 
and  'twist  design'  respectively.  In  the  stress  design,  the  stresses  in 
all  the  elements  are  less  than  or  equal  to  the  maximum  allowable  stress 
and  in  the  twist  design,  the  twist  of  the  tip  of  the  wing  is  equal  to 
the  prescribed  value. 
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The  iteration  history  for  the  negative  twist  constraint  problem  is 

given  in  Table  8.  For  the  first  three  iterations  A . < A , hence  the 

d s 

weights  corresponding  to  the  stress  design  are  larger  than  the  twist 
design.  The  weight  of  the  structure  is  directly  proportional  to  the 
scaling  parameter.  The  absolute  value  of  the  twist  of  the  tip  of  the 
wing  for  the  stress  design  is  less  than  7°  for  the  first  three  iterations. 
After  the  fourth  iteration  > As  and  thus  the  weight  of  the  twist 
design  is  greater  than  the  stress  design.  For  these  iterations  the  twist 
design  has  the  tip  rotation  of  -7°  and  the  stresses  in  the  elements  are 
less  than  the  maximum  allowables.  The  least  weight  design  of  71.86  lbs  is 
obtained  at  the  11th  iteration.  The  weights  for  the  stress  design  and 
the  twist  design  are  plotted  for  each  iteration  in  Figure  6.  This 
figure  shows  that  after  the  fourth  iteration,  there  is  no  significant 
difference  in  the  weight  of  the  structure  satisfying  the  twist  constraint. 
The  design  has  reached  a stable  position  in  the  design  space.  The 
displacements  at  nodes  7 and  9 for  the  least  weight  design  are  given  in 
Table  9.  The  thicknesses  and  cross-sectional  areas  of  members  for  the 
minimum  weight  design  are  given  in  Tables  10  and  11.  The  distribution 
of  the  number  of  plies  in  the  four  fiber  directions  for  the  minimum 
weight  design  is  given  in  Figure  7. 

For  the  positive  twist  constraint  of  +2°  (wash  in),  the  iteration 
history  is  given  in  Table  12.  The  twist  of  the  tip  for  the  initial 
design  is  negative,  and  this  design  cannot  be  scaled  to  achieve  +2" 
twist.  From  the  2nd  to  the  8th  iteration  the  twist  of  the  stress  design 
increases  with  each  iteration  until  it  reaches  +2°.  In  subsequent 
iterations  for  nearly  ail  cases,  the  twist  designs  are  acceptable  designs 
since  A^  > A . The  weight  of  the  structure  versus  the  iteration  numbers 
is  plotted  in  Figure  8.  Figure  8 shows  that  the  weight  of  the  structure 
oscillates  and  the  design  does  not  reach  a stable  position  in  the  design 
space.  For  the  twenty-one  iterations,  the  least  weight  design  is 
obtained  in  the  last  iteration.  The  displacements  at  nodes  7 and  9 for 
this  design  are  given  in  Table  13.  Tables  14  and  15  contain  the 
details  of  the  design.  The  number  of  plies  in  the  top  and  bottom  skin 
is  given  in  Figure  9. 

13 
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SECTION  V 
CONCLUSION 


It  has  been  demonstrated  that  an  algorithm  based  on  optimality 
criteria  can  be  used  to  design  a minimum  weight  wing  structure  with 
specific  negative  (wash  out)  or  positive  (wash  in)  twist.  The  design 
satisfies  the  twist  constraint  and  specific  strength  criteria  in  each 
element. 

Future  work  in  this  area  will  be  devoted  to  developing  a method 
to  design  a wing  type  structure  with  specified  twist  at  more  than  one 
location.  Further  work  is  also  needed  to  develop  a program  to  design 
a wing  structure  for  twist,  strength,  and  stiffness  for  flutter  and 
divergence  requirements. 


Figure  3.  Type  of  Elements 


30  3° 


Model  of  Wing 


NUMBER  OF  LAYERS  WITH  FIBER  ORIENTATIONS  IN  0%  90°,  +45°,  -45°  DIRECTION  RESPECTIVELY 


Design  of  Wing  (Stress  Constra 


Figure  6.  Iteration  for  Twist  Constraint  (Wash  out 
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Figure  7.  Optimum  Design  of  Wing  (Twist  Constraint,  Wash  out 


Stress  Design 
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TABLE  1 

ELEMENT  CONNECTIONS* 


Element 

No. 

A 

B 

C 

D 

Element 

No. 

A 

B 

C 

D 

1 

1 

9 

11 

C 

61 

73 

73 

85 

83 

2 

2 

13 

12 

7 

62 

74 

7 6 

86 

84 

3 

1 

3 

1 3 

1 1 

63 

75 

77 

8 7 

85 

4 

2 

4 

1 4 

12 

6 4 

76 

78 

9 « 

4b 

5 

3 

5 

1 5 

13 

6 5 

1 

2 

1 0 

c 

6 

4 

6 

1 6 

1 4 

66 

1 

2 

4 

7 

7 

5 

7 

17 

15 

67 

3 

4 

6 

c 

9 

6 

9 

1 3 

16 

69 

5 

6 

8 

7 

9 

9 

11 

21 

19 

69 

9 

10 

1 2 

1 1 

n 

13 

12 

22 

20 

7 : 

11 

12 

14 

13 

11 

11 

13 

2 3 

21 

71 

13 

Is 

1 6 

lc 

12 

12 

14 

24 

22 

72 

15 

16 

1 4 

1 7 

13 

13 

15 

2 5 

23 

73 

19 

23 

22 

21 

14 

14 

1 o 

26 

24 

74 

21 

2 7 

2s 

23 

15 

15 

17 

27 

25 

75 

23 

24 

2 

25 

16 

16 

19 

23 

26 

76 

25 

26 

2 9 

27 

17 

19 

21 

31 

29 

77 

29 

3 0 

39 

31 

18 

20 

22 

32 

3 0 

79 

31 

32 

3 s 

33 

19 

21 

2 3 

3 3 

31 

79 

33 

34 

35 

35 

2: 

22 

24 

3 4 

3 2 

80 

35 

3 6 

3 8 

37 

21 

23 

25 

3 5 

3 3 

81 

39 

s3 

4 -> 

si 

22 

24 

2 6 

36 

34 

82 

41 

42 

4 s 

s3 

23 

25 

27 

37 

35 

8 3 

43 

4 4 

46 

45 

24 

26 

29 

3 9 

76 

84 

45 

4 6 

44 

47 

25 

29 

31 

41 

39 

8 5 

49 

53 

5 7 

C1 

26 

33 

32 

42 

4 C 

86 

51 

52 

t?  4 

53 

27 

31 

33 

4 3 

41 

87 

53 

54 

5 6 

c,  r. 

28 

32 

34 

44 

*2 

88 

c 5 

56 

5 4 

5 7 

29 

33 

35 

4 5 

43 

89 

59 

60 

6? 

61 

3: 

34 

36 

4 6 

s4 

9 

61 

62 

6 ^ 

6 3 

31 

35 

37 

4 7 

45 

91 

6 3 

6 4 

b 0 

65 

32 

36 

38 

4 9 

46 

92 

6 5 

66 

64 

67 

33 

39 

41 

51 

49 

93 

69 

7 0 

7 2 

71 

34 

40 

42 

52 

50 

94 

71 

72 

7 -♦ 

7 7 

35 

41 

43 

5 3 

51 

9 5 

73 

74 

76 

75 

36 

42 

4*4 

54 

52 

96 

75 

7 3 

7 * 

77 

37 

43 

4-3 

55 

u 3 

97 

9 

10 

2 3 

19 

39 

44 

46 

56 

6 4 

OH 

19 

20 

30 

29 

39 

45 

4 7 

^ 7 

5 6 

99 

29 

i 3 

4 7 

3 9 

4 0 

46 

48 

59 

5 6 

10  3 

39 

4 0 

50 

- 9 

41 

49 

51 

61 

59 

101 

49 

53 

6 1 

59 

42 

5 3 

52 

6 2 

60 

102 

59 

60 

’I 

6 c 

43 

51 

53 

63 

61 

ir  3 

69 

7 0 

8 T 

79 

44 

52 

54 

6 4 

62 

10  4 

3 

4 

1 s 

1 3 

45 

53 

55 

65 

63 

1 0 5 

13 

Is 

24 

2 3 

46 

54 

56 

6 6 

6 4 

10  6 

2 3 

24 

34 

3 7 

47 

55 

57 

67 

65 

1 C 7 

93 

34 

4 -4 

43 

4 R 

56 

58 

68 

6 6 

108 

43 

4 4 

5 4 

53 

49 

59 

61 

7 1 

69 

10  9 

53 

5s 

b 4 

6 3 

50 

63 

62 

72 

70 

11  " 

6 3 

6 4 

74 

73 

51 

61 

6 3 

7 3 

71 

1 1 1 

73 

7s 

4 4 

6 3 

52 

62 

6 4 

74 

72 

112 

7 

9 

1 4 

1 7 

53 

63 

65 

75 

73 

113 

17 

1 9 

2 4 

2 7 

54 

64 

66 

76 

74 

114 

27 

24 

3 4 

3 7 

55 

65 

67 

77 

75 

115 

3 7 

38 

8 

- 7 

56 

66 

6 3 

7 9 

76 

1 1 A 

47 

44 

4 < 

6 7 

57 

69 

71 

9 1 

79 

11  7 

5 7 

54 

6 « 

6 7 

5 9 

73 

72 

8 2 

80 

118 

67 

6 4 

7 4 

7 7 

59 

71 

73 

4 4 

4 1 

119 

77 

7 9 

4 4 

6 7 

60 



72 

74 

8s 

82 





*See  Figures  3 and  4 


Elements  120-158  are  posts  connecting  the  top  and  bottom  nodes. 
For  example,  element  120  connects  nodes  1 and  2 and  element  158 
connects  nodes  77  and  78. 
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TABLE  2 (Concluded) 


FO^CE- Y 


39J.  I 
4 0 *4  « • 
<.0t.  I 

423  . : 
420.  i 
1120. 
1120.  ! 
6 8 3. 

8 8 3.1 
412. 

413.1 
391 . 
391. 
36*.  i 
36  *.  i 
6 0'..! 
8t  4.  : 

1 0 40  . 
1C40. 
433. 

4 3 3. 

370.1 
27  0 . 

304.1 

3 u 

4 4 6 

4 4 6 

5 4 1. 
64-1 
224 
224 
190 
190 
159 

1 55. 
??(- 
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TABLE  3 

PROPERTIES  AND  ALLOWABLE  STRESSES 


! MATERIAL  ! 

Properties 

Graphite  Epoxy 

A1 uminum 

En  (psi) 

18.5  x 106 

10.5  x 10fc 

E22  (psi) 

1 .6  x 106 

10.5  x 106 

Poisson's  Ratio 

0.25 

0.3 

Shear  Modulus  (psi ) 

0.65  x 106 

4.038  x 106 

Specific  Weight 

0.055 

0.1 

(lbs/in3) 

Lamina  Thickness( in) 

0.0052 

— 

Allowable  Stresses  (Ksi) 

F (tension) 

139.0 

45.0 

Fx  (compression) 

86.0 

45.0 

F (tension) 

y 

— 

45.0 

Fy  (compression) 

— 

45.0 

F 

46.8 

25.9 

xy 
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TABLE  5 

ELEMENT  AND  LAYER  THICKNESSES  FOR  THE 
MINIMI 'M  WEIGHT  STRESS  CONSTRAINT  DESIGN 


Thickness  (in.) 

Element  No. 

Total 

0° 

O 

O 

O' 



. 

■O 

o 

+45° 

1 

.0418 

.0104 

.0104 

.0106 

.0104 

2 

.0418 

.0104 

.0104 

.0106 

.0104 

3 

.0484 

.0108 

.0108 

.0159 

.0108 

4 

.0420 

.0104 

.0104 

.0108 

.0104 

5 

.0419 

.0107 

.0104 

.0104 

.0104 

6 

.0419 

.0107 

.0104 

.0104 

.0104 

7 

.0643 

.0189 

.0104 

.0246 

.0104 

8 

.0477 

.0117 

.0104 

.0152 

.0104 

9 

.0420 

.0108 

.0104 

.0104 

.0104 

10 

.0470 

.0103 

.0104 

.0104 

.0104 

11 

.0430 

.0118 

.0104 

.0104 

.0104 

12 

.0419 

.0107 

.0104 

.0104 

.0104 

13 

.0563 

.0251 

.0104 

.0104 

.0104 

14 

.0425 

.0113 

.0104 

.0104 

.0104 

15 

.0548 

.0216 

.0104 

.0124 

.0104 

16 

.0426 

.0114 

.0104 

.0104 

.0104 

17 

.0423 

.0111 

.0104 

.0104 

.0104 

18 

.0421 

.0109 

.0104 

.0104 

.0104 

19 

.0723 

.0387 

.0112 

.0112 

.0112 

20 

.0465 

.0153 

.0104 

.0104 

.0104 

21 

.0697 

.0385 

.0104 

.0104 

.0104 

22 

.0507 

.0195 

.0104 

.0104 

.0104 

23 

.0765 

.0419 

.0106 

.0135 

.0106 

24 

.0576 

.0247 

.0110 

.0110 

.0110 

25 

.0426 

.0114 

.0104 

.0104 

.0104 

26 

.0422 

.0110 

.0104 

.0104 

.0104 

27 

.0974 

.0628 

.0115 

.0115 

.0115 

28 

.0598 

.0285 

.0104 

.0104 

.0104 

29 

.0971 

.0659 

.0104 

.0104 

.0104 

30 

.0641 

.0329 

.0104 

.0104 

.0104 

31 

.0908 

. 0566 

.0104 

.0135 

.0104 

32 

.0719 

.0386 

.0111 

.0111 

.0111 

33 

.0426 

.0114 

.0104 

.0 104 

.0104 

34 

.0422 

.0110 

.0104 

.0104 

.0104 

35 

.1180 

.0831 

.0116 

.0116 

.0116 

36 

.0725 

.0409 

.0105 

.0105 

.0105 

37 

.1320 

.0978 

.0106 

.0106 

.0130 

38 

.0806 

.0494 

.0104 

.0104 

.0104 

39 

. 1040 

.0716 

.0104 

.0117 

.0104 

40 

.0829 

.0508 

.0107 

.0107 

.0107 

41 

.0425 

.0113 

.0104 

.0104 

.0104 

42 

.0422 

.0110 

.0104 

.0104 

.0104 

43 

. 1 340 

.0987 

.0118 

.0118 

.0118 

44 

.0822 

.0503 

.0106 

.0106 

.0106 

45 

.1710 

.1260 

. 0 106 

.0106 

.0236 

2 R 


TABLE  5 (Concluded) 


Thickness  (in. 

) 

Element  No. 

Total 

0° 

90° 

-45° 

+45° 

46 

.0999 

.0665 

.0104 

.0104 

.0126 

47 

.1390 

.1078 

.0104 

.0104 

.0104 

48 

. 1010 

.0697 

.0104 

.0104 

.0104 

49 

.0424 

.0112 

.0104 

.0104 

.0104 

50 

.0422 

.0110 

.0104 

.0104 

.0104 

51 

.1360 

.1017 

.0115 

.0115 

.0115 

52 

.0855 

.0537 

.0106 

.0106 

.0106 

53 

.2100 

.1431 

.0107 

.0107 

.0456 

54 

.1230 

.0774 

.0104 

.0104 

.0249 

55 

. 1810 

.1498 

.0104 

.0104 

.0104 

56 

.1270 

.0958 

.0104 

.0104 

.0104 

57 

.0421 

.0109 

.0104 

.0104 

.0104 

58 

.0421 

.0109 

.0104 

.0104 

.0104 

59 

.1580 

.1192 

,0109 

.0109 

.0170 

60 

.0969 

.0655 

.0105 

.0105 

.0105 

61 

.2300 

.1645 

.0107 

.0107 

.0441 

62 

.1360 

.0907 

.0104 

.0104 

.0245 

63 

.2510 

.1828 

.0104 

.0104 

.0474 

64 

.1690 

.1175 

.0104 

.0104 

.0307 
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TABLE  6 

SHEAR  PANEL  AND  BAR  DIMENSIONS  FOR  THE 
MINIMUM  WEIGHT  STRESS  CONSTRAINT  DESIGN 


Shear  Panels 


Element  No. 


65-104 

105 

106 

107 

108 

109 

110 
111 


Thickness  (in.) 


.0416 

.0446 

.0658 

.0827 

.0900 

.0926 

.0992 

.1090 


Shear 

Panels 

Element  No. 

Thickness  (in.) 

112 

.0416 

113 

.0515 

114 

.0570 

115 

.0637 

116 

.0772 

117 

.0973 

118 

.1210 

119 

.1390 

Bars 

Element  No. 

Area  (in^) 

120-158 

.0416 

TABLE  7 

DEFLECTIONS  AT  TIP  FOR  THE 
MINIMUM  WEIGHT  STRESS  CONSTRAINT  DESIGN 


Node  7 (trailing  edge) 
Direction  (in) 

Node  9 (leading  edge) 
Direction  (in) 

Twist  (deg.) 

X 

Y 

Z 

X 

Y 

Z 

-.0565 

-. 1961 

13.4306 

-.0435 

-.1632 

10.6296 

-5.458 

30 
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TABLE  8 

ITERATIVE  HISTORY  FOR  THE  TWIST  CONSTRAINT  PROBLEM  (WASH  OUT  7#) 


Cycle 

Weight 

(lbs.) 

Twist 

(deg. ) 

No. 

M-n  mg 

■ssssmssH 

Stress  Design 

Twist  Design 

— 

309.34 

82.70 

-1.879 

-7.000 

213.96 

153.70 

-5.040 

-7.000 

76.61 

46.31 

-4.244 

-7.000 

58.72 

72.22 

-8.589 

-7.000 

63.49 

73.87 

-8.130 

-7.000 

63.89 

76.83 

-8.400 

-7.000 

63.78 

72.70 

-7.967 

-7.000 

8 

63.71 

76.82 

-8.421 

-7.000 

HOB 

62.95 

72.93 

-8.097 

-7.000 

■1 

64.56 

76.51 

-8.280 

-7.000 

64.28 

71.86 

-7.815 

-7.000 

63.40 

76.07 

-8.381 

-7.000 

63.53 

72.35 

-7.959 

-7.000 

! 

I 


TABLE  9 

DEFLECTIONS  AT  TIP  FOR  THE  MINIMUM  WEIGHT 
TWIST  CONSTRAINT  DESIGN  (WASH  OUT  7*) 


Node  7 (trailing  edge) 
Direction  (in) 

Node  9 (leading  edge) 
Direction  (in) 

Twist  (deg.) 

X 

Y 

Z 

X 

Y 

Z 

-.0995 

-.1594 

11.6800 

-.0881 

-.1073 

8.0799 

-7.00 

31 
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TABLE  10 


ELEMENT 


AND  LAYER  THICKNESSES  FOR  THE  MINIMUM  WEIGHT 
TWIST  CONSTRAINT  DESIGN  (WASH  OUT  7°) 


Thickness  (in. 

) 

Element  No. 

Total 

0° 

90° 

-45° 

+45° 

1 

.0494 

.0124 

.0124 

.0124 

.0124 

2 

.0494 

.0124 

.0124 

.0124 

.0124 

3 

.05  79 

.0208 

.0124 

.0124 

.0124 

4 

.0494 

.0124 

.0124 

.0124 

.0124 

5 

.0494 

.0124 

.0124 

.0124 

.0124 

6 

.0494 

.0124 

.0124 

.0124 

.0124 

7 

.0736 

.0191 

.0123 

.0298 

.0123 

8 

.0554 

.0123 

.0123 

.0184 

.0123 

9 

.1550 

.0123 

.0302 

.0123 

.1001 

10 

.1360 

.0124 

.0250 

.0124 

.0863 

11 

.0572 

.0201 

.0124 

.0124 

.0124 

12 

.0494 

.0124 

.0124 

.0124 

.0124 

13 

.'.678 

.0307 

.0124 

.0124 

.0124 

14 

.0494 

.0124 

.0124 

.0124 

.0124 

15 

.0530 

.0153 

.0124 

.0130 

.0124 

16 

.0494 

.0124 

.0124 

.0124 

,0124 

17 

.1940 

.0124 

.0513 

.0124 

.1179 

18 

.1250 

.0124 

.0338 

.0124 

.0665 

19 

.0988 

.0617 

.0124 

.0124 

.0124 

20 

.0590 

.0220 

.0123 

.0123 

.0123 

21 

.0768 

.0397 

.0124 

.0124 

.0124 

22 

.0581 

.0210 

.0124 

.0124 

.0124 

23 

.0769 

.0347 

.0124 

.0175 

.0124 

24 

.0595 

.0225 

.0123 

.0123 

.0123 

25 

.2020 

.0124 

.0538 

.0124 

.1234 

26 

.1330 

.0124 

.0341 

.0124 

.0742 

27 

.1400 

.0818 

.0123 

.0123 

.0335 

28 

.0937 

.0298 

.0124 

.0124 

.0392 

29 

.1240 

.0581 

.0124 

.0124 

.0411 

30 

.0947 

.0290 

.0124 

.0124 

.0410 

31 

.0962 

.0591 

.0124 

.0124 

.0124 

32 

.0769 

.0398 

.0124 

.0124 

.0124 

33 

.2290 

.0123 

.0699 

.0124 

.1344 

34 

.2190 

.0124 

.0719 

.0124 

.1224 

35 

.1750 

.0947 

.0123 

.0123 

.0551 

36 

. 1110 

.0381 

.0124 

.0123 

.0482 

37 

. 1560 

.0610 

.0123 

.0123 

.0703 

38 

.1020 

.0309 

.0124 

.0124 

.0464 

39 

.1500 

.1129 

.0124 

.0124 

.0124 

40 

. 1090 

.0719 

.0124 

.0124 

.0124 

41 

.1580 

.0124 

.0124 

.0124 

.1209 

42 

.1830 

.0124 

.0209 

.0124 

.1373 

43 

.2120 

.0823 

.0123 

.0123 

.1051 

44 

.1560 

.0289 

.0124 

.0124 

.1024 

45 

. 1960 

.0536 

.0123 

.0123 

.1177 

3? 
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TABLE  10  (Concluded) 


Element  No. 


Total 


Thickness  (in.) 


.1480 
.2250 
.1810 
.0678 
.0780 
.2460 
.2020 
.2830 
.2300 
.3200 
.2270 
.0780 
.0780 
.0759 
.0704 
.4040 
.3800 
. 3970 
.2700 


.0123 

.0124 

.0124 

.0123 

.0123 

.0123 

.0123 

.0124 

.0124 

.0124 

.0124 

.0123 

.0123 

.0124 

.0123 

.0124 

.0124 

.0123 

.0123 


024ft  | 

0410 
1714 
1650 
2162 
19  30 

0410 

0411  | 
0123 
0123 
0388 

0 3 34 
3228 
3430 
0472 
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TABLE  11 

SHEAR  PANEL  AND  BAR  DIMENSIONS  FOR  THE  MINIMUM  WEIGHT 
TWIST  CONSTRAINT  DESIGN  (WASH  OUT  7*) 


Shear  Panels 

Element  No.  Thickness  (in.) 


Shear 

Panels 

Element  No. 

Thickness  (in.) 

96 

.0485 

97 

.1380 

98 

.1370 

99 

.1490 

100 

.1450 

101 

.0898 

102-105 

.0485 

106 

.0605 

107 

.0552 

108 

.0485 

109 

.0655 

110 

.1340 

111 

.1950 

112 

.0485 

113 

.0592 

114 

.0681 

115 

.0968 

116 

.1340 

117 

.1690 

118 

.2070 

119 

.2290 

Node  No. 


Area  (in^) 


Node  No. 


Area  (in^) 


20-122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

33-137 

38-140 

141 


142 

143 

144 

145-147 

148 

149 

150 

151 

152 

153 

154 

155 

156-158 


34 
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TABLE  12 

ITERATIVE  HISTORY  FOR  THE  TWIST  CONSTRAINT  PROBLEM  (WASH  IN  2°) 


Cycle 

Weight 

(lbs.) 

Twist 

(deg. ) 

No. 

Stress  Design 

Twist  Design 

Stress  Design 

Twist  Design 

1 

309.34 

-1.879 

2 

434.58 

8.24 

.0379 

2.000 

3 

558.98 

116.80 

.418 

2.000 

4 

461.96 

122.51 

.531 

2.000 

5 

422.68 

230.15 

1.090 

2.000 

6 

546.74 

515.48 

1.886 

2.000 

7 

370.82 

253.58 

1. 368 

2.000 

8 

416.46 

416.46 

2.000 

2.000 

9 

258.97 

258.97 

2.000 

2.000 

10 

326. 38 

326.38 

2.000 

2.000 

11 

135.14 

126.67 

1.875 

2.000 

12 

164.13 

164.13 

2.000 

2.000 

13 

194.59 

194.59 

2.000 

2.000 

14 

230.57 

230.57 

2.000 

2 . 000 

15 

206.63 

206.63 

2 . 000 

2.000 

16 

248.42 

248.42 

2.000 

2.000 

17 

127.60 

127.60 

2.000 

2.000 

18 

274.28 

274.28 

2.000 

2.000 

19 

232.19 

232.19 

2.000 

2.000 

20 

260.97 

260.97 

2.000 

2.000 

21 

105.54 

105.54 

2.000 

2.000 

TABLE  13 

DEFLECTIONS  AT  TIP  FOR  THE  MINIMUM  WEIGHT 
TWIST  CONSTRAINT  DESIGN  (WASH  IN  2°) 


Node  7 (trailing  edge) 
Direction  (In) 

Node  9 (leading  edge) 
Direction  (in) 

Twist  (deg.) 

X 

Y 

z 1 

X 

Y 

Z 

.0607 

-.1073 

4 . 49  31 

.0449 

-.1172 

5.5179 

2 . 000 

35 
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TABLE  14 

ELEMENT  AN1)  LAYER  TH1CKNESSESS  FOR  THE  MIN IMl'M  WI  ICliT 
TWIST  CONSTRAINT  DESIGN  (WASH  IN  2°) 


Thickness  (in. 

2 ' . 1 

Element  No. 

total 

0° 

90° 

-45° 

+45° 

1 

.0680 

.0170 

.0170 

.0170 

.0170 

2 

.0680 

.0170 

.0170 

.0170 

.0170 

3 

.1290 

.0778 

.017i 

.0171 

.0171 

4 

.1280 

.0769 

.0170 

.0170 

.0170 

5 

.0680 

.0170 

.0170 

.0170 

.0170 

6 

.0708 

.0170 

. 0 1 70 

.0197 

.0170 

7 

. 3880 

.2877 

.0170 

.06  38 

.0194 

8 

. 3870 

. 2807 

.0170 

.0723 

.0170 

9 

.0680 

.0170 

.0170 

.0170 

.0170 

10 

.0680 

.0170 

.0170 

.0170 

.0170 

11 

.0680 

.0170 

.0170 

.0170 

.0170 

12 

.0680 

.0170 

.0170 

.0170 

.0170 

13 

.0695 

.0170 

.01  70 

.0170 

.0185 

14 

. 0680 

.01  70 

.0170 

.0170 

.0170 

15 

.5460 

.4290 

.0170 

.0829 

.0170 

lb 

.5400 

.425  3 

.0170 

.0806 

.0170 

17 

.0680 

.0170 

.0170 

.0170 

.0170 

18 

.0680 

.0170 

.0170 

.0170 

.0170 

19 

.0680 

.01  70 

.0170 

.01  70 

.0170 

20 

.0680 

.01  70 

.01 

.(! : 70 

.0170 

) 1 

.0695 

.0170 

.0170 

.0170 

.0185 

22 

.0680 

. 0 1 70 

.0170 

.0170 

. 170 

23 

. 7190 

.5  308 

.0611 

. 1 102 

.0170 

24 

.6830 

.5167 

.0352 

.1141 

.0170 

25 

.0680 

. 01  70 

.0170 

.0170 

.0170 

26 

.0680 

.0170 

.0170 

. 0 i 70 

.0170 

2 7 

.0680 

.01  70 

.0170 

.0170 

.0170 

28 

.0680 

.0170 

.0170 

.0170 

. 

1 29 

.1420 

.0905 

.0170 

.iii  70 

.0174 

30 

.0828 

.0318 

.0170 

.0170 

.0170 

31 

.6940 

.4272 

.0561 

. 1936 

.0170 

32 

. 7660 

.4265 

. 1092 

.2133 

.0170 

33 

.0680 

.0170 

.0170 

.0170 

.0170 

34 

.0680 

.01 70 

.0170 

.0170 

.0170 

35 

.0680 

.0170 

.0170 

.0170 

.0170 

36 

.0680 

.01  70 

.0170 

.0170 

.0170 

37 

.0850 

.0170 

. 1170 

.0170 

.0339 

38 

.0692 

.0170 

. 03  70 

. 0 1 82 

.01  70 

39 

. 7290 

.6780 

. 0 1 70 

.01  70 

.0170 

40 

.6900 

.6390 

.01  70 

. 0 1 70 

.0170 

41 

. 0680 

.0170 

.01/0 

.0170 

.0170 

1 42 

.0680 

.0170 

.0170 

.0170 

.0170 

43 

.0680 

.01  70 

.0170 

.0!  70 

.0170 

44 

. ;m.80 

.0170 

.0170 

.01  7') 

.0170 

1 45 

i 

.0919 

i 

. M 70 

.01  70 

.01  70 

.0408 
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TABLE  14  (Concluded) 


Thi.'kiiess  (in.) 

Element  No. 

Total 

. ___  o°  . 

90° 

-45° 

445° 

46 

.0767 

.0170 

.0170 

.0256 

.0170 

47 

.7790 

. 7280 

.0170 

.0170 

.0170 

48 

.7400 

.6889 

.0170 

.0170 

.0170 

49 

.0680 

.0170 

.iM  70 

.01  70 

. 01 70 

50 

.0680 

.0170 

.0170 

.0170 

.0170 

51 

.0709 

.0199 

.0170 

.0170 

.0170 

52 

.0680 

.0170 

.0170 

.0170 

.0170 

53 

.0927 

.0170 

.0170 

.0170 

.0416 

54 

.7290 

. 1701 

.1701 

. J701 

.2188 

55 

.7870 

. 7359 

.0170 

.0170 

.0170 

56 

.7780 

.7270 

.0170 

.0170 

.0170 

57 

.0680 

.0170 

.0170 

.0170 

.0170 

58 

.0680 

.0170 

.0170 

.0170 

.0170 

59 

.0807 

.0297 

.01  70 

.0170 

.0170 

60 

. 0680 

.0170 

.0170 

.0170 

.0170 

61 

. 1 300 

.0755 

.0170 

.0170 

.0206 

62 

.0852 

.0342 

.01  70 

.0170 

.0170 

63 

.7650 

.7140 

.0170 

.0170 

.0170 

64  

.7800 

. . ■ 

.0170 

.0170 

.0170 

1 
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